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ABSTRACT 

We report the detection of Galactic supernova remnants (SNRs) in the mid-infrared (at 24 and 70 
/im), in the coordinate ranges 10° < I < 65° and 285° < / < 350°, \h\ < 1°, using the Multiband 
Imaging Photometer (MIPS) aboard the Spitzer Space Telescope. We search for infrared counterparts 
to SNRs in Green's catalog and identify 39 out of 121, i.e., a detection rate of about 32%. Such a 
relatively low detection fraction is mainly due to confusion with nearby foreground/background sources 
and diffuse emission. The SNRs in our sample show a linear trend in [F8/i^24] versus [F70/-F24]. We 
compare their infrared fluxes with their corresponding radio flux at 1.4 GHz and find that most 
remnants have ratios of 70 /im to 1.4 GHz characteristic of SNRs (with the exception of a few which 
have ratios closer to those of H II regions). Furthermore, we retrieve a slope close to unity when 
correlating infrared (24 and 70 /im) with 1.4 GHz emission. Our survey is more successful in detecting 
remnants with bright X-ray emission, which we find is well correlated with the 24 /Ltm morphology. 
Moreover, by comparing the power emitted in X-ray, infrared and radio, we conclude that the energy 
released in the infrared is comparable to the cooling in the X-ray range. 
Subject headings: infrared: ISM — shock waves — supernova remnants 



1. INTRODUCTION 



With the h elp of Galactic surveys (e.g., Arendt 1989 
Saken et al. 1992^), our understanding of the inirarec 



(IR) energetics, morphology and evolution of supernova 
remnants (SNRs) has increased substantially in the past 
decades. However, even at IRAS (Infrared Astronomical 
Satellite) resolution, confusion with other IR sources in 
the Galactic plane is a limiting factor that can prevent 
a clear picture from emerging. Such confusion, for ex- 
ample, can make the correct assessment of ejecta dust 
masses for Type II SNRs difficult, a key factor for in- 
terpreting high redshift dust in the Universe (ISibthorpej 
et al.|i!2010, and references therein). Infrared emission 
from SNRs can also provide insight i nto the radio -IR cor- 
relation of external galaxies (Hclou et al. 1985). Addi- 
tionally, the study of dust re-emission is relevant for de- 
termining the cooling rate of SN Rs, an impo rtant quan- 
tity that shapes their evolution ( Dwek|[T987 ) . 

In this paper, we search for counterparts of SNRs in 
the mid-infrared (at 24 and 70 /im) using the higher 
resolution Spitzer data obtained with the Multiband 
Imaging Photometer (MIPS), in the coordinate ranges 
10° < ; < 65° and 285° < / < 350°, I6| < 1°, comple- 
menting the work of iReach et al. (2006). We also briefly 



explore the relation of SNRs IR fluxes with other wave- 
length ranges, such as the radio and the X-ray. 

The structure of this paper is as follows: In §1.1[ we 
review previous Galactic infrared surveys and their de- 
tection rates. In p.2[ we discuss SNR emission mecha- 
nisms which can contribute in the mid-infrared regime. 
In 32] we describe the mid-infrared, radio and X-ray data 
usea in this work. In fj3] we report the MIPS detections. 
A discussion of SNR identifications, derived color ratios 
and radio and X-ray to infrared ratios follows in ^ Fi- 
nally, in ^ we summarize our findings. Details on each 



detected SNR can be found in the appendix. 

1.1. Previous Galactic Infrared Surveys 

Using the Galactic Legacy In frared Mid-Plane Sur 
vey Extraordinaire, GLIMPSE (jBeniamin et al.||2003 



ChurchweU et al.|2009[ ), [Reach et al.| ( |2006[ ) searcKe^TEe 



Galactic Plane for mtrared counterparts to known SNRs. 
Out of 95 objects, 18 were clearly detected, a detection 
rate of about 20%. Previous searches were conducted 
using the IRA S all-sky surve y. Within the GLIMPSE 



surveyed ar ea, Arendt 



of 17% and Saken et a. 



( 1989 ) obtain ed a detection rat e 
| ( |iy92| ) 18% dReach et al.||2006 ). 
In the present study, we identify 39 out of 121, i.e., a 
detection rate of about 32%. 

1.2. SNR infrared emission mechanisms 

The IR part of the spectral energy distribution (SED) 
of SNRs can result from multiple sources ranging 
from dust (either stochastically or thermally heated), 
atomic/molecular line emission, PAH (polycyclic aro- 
matic hydrocarbon) emission, and even synchrotron 
emission. In fact, the amount of the latter emission 
can vary substantially and is dependent on the type of 
remnant. In plerion rem nants such as the Crab Nebula 
( Strom fc Greidanus|1992| ), it accounts for about 90% of 
the mid-mtrared flux density (at 25 /um). In contrast, 
about 2% of the mid-infrared flux density (at around 25 
lim) in Cas A is due to synchrotron emission JSibthorpe 
eral][2010| . 



Dust emission is expected to be a substantial com- 
ponent of the IR emission of SNRs since their spectra 
is well fitted by one or more thermal dust pop ulations, 
as revealed by previous infrared surveys (IRAS; Arendt 
1989 Saken et al.||1992), or by stoc hastic heating of the 
grams, as proposed by Hines et al. ( 2004j ) to explain the 



2 



mid-infrared emission from Cas A. The dust is heated 
by c harged particles in hot plasma generated by shocks 



some Galactic SNRs has also been shown throug h near 



e^ |Dwek fc Werner|198l||Dwekll987[|Dwek fc Arendt 
9921). Direct evidence tor the snock-neated plasma is 



1992 1^ Direct evidence tor the snock-neated plasma is 
provided by X-ray observations in the continuum, and 
X-ray lines are a good diagnostic o f high energy inter- 
actions and abundances ( Vink||2004 ) . While for the ma- 
jority of SNRs, the observed dust emission is thought 
to originate from the interaction of the Shock wave with 
the su rrounding interstellar mediu m (ISM; e.g., Douvion 
erar][20 01, Bor kowski ct al. 20061 | Williams et aLpOOnaT 
Blair et al. || 2 007 j, evidence of SNRs with ejecta dust is 
scarce (e.g., Arendt et al.| 1999) even with the help of 
higher resolution IR telescopes such as ISO and Spitzer. 
However, obs ervations with t he Spitzer Infrared Spectro- 
graph (IRS) ( Rho et al.|2008 ) of the young remnant Cas 
A have pointed to the presence of a large hot dust contin- 
uum peaking at 21 /xm, from ejecta dust formed through 
condensation of freshly synthesized heavy elements. This 
component has temperatures ranging from 60 to 120 K. 
Furthermore, recent far-infrared and sub- millimeter work 



(Sibthorpe et al. 20101 on the same remnant, also sug- 
gested the existence of a 'tepid' 35 K) central dust 
comp onent in the ejected material (see also Barlow et al. 

Mo| . 

Observational evidence for the importance of IR lines 
has been provided by multiple SNRs spectroscopic stud- 
ies which identifie d a plethora of lines from elements such 
as Fe and O (e.g., |01iva et a l."1999bl |Reach et al]|2002[ ). 
Much of the infrared line emission seen in SIN Rs results 
from the interaction of the blastwave with the surround- 
ing environment where they are born as the shocked 
gas cools down. For example, IC443 is encountering an 
atomic region at the northeas t side and a molecular cloud 
at the southern border (e.g. ,|Rho et al. 2001 ). Spectral 
observations by Oliva et al. (1999a) on the north of the 
remnant revealed that most of the emission at 12 and 25 
/xm, previously surveyed with IRAS, was in fact due to 
ionized line emission (e.g., [Fe II]) with only a small con- 
tribution from dust. Because optical spectra also showed 
strong coUisionally excited fine structure atomic lines 
( [Fesen fc Kirshner^l980j , the infrared line emission is not 
unexpected. Likewise, for the 'optically bright regions' 
in N49, an old remnant in the Large Magellanic Cloud 
(LMC), [wriliams et al |^2006b[ ) found that IR line emis- 
sion can be a substantial traction of the total IR flux, up 
to 80% at 24 ^m. On the south of IC443, however, the 
spectrum is domi nated by the Ho pur e rotational lines 
S(2) through S(7) ( Cesarsky et al.|l99 9 [), once aga i n with 
a faint dust continuum. Moreover, Troja et al. (20061 
found a good agreement between the soft X-ray emitting 
plasma and the radio/optical structure on the northeast- 
ern part of the remnant where the plasma density is high- 
est. Such energetic encounters change the morphology of 
the neighboring ISM and of the remnant itself thereby 
creating prominent emission lines (e.g., molecular hydro- 
gen) which can be detected in the infrared. OH masers 
(in particular, the 1720 MHz line) are also good tracers 
of the remnant encounter with dense molecular regions 
with 10% of Galactic S NRs showing associated masers 



( Wardle &: Yusef-Zadehj|2002| and references therein). 
I'he existence ol infrared emission associated wi 



with 

PAHs in SNR s was first observed in a LMC remnant by 
Tappe et al. ( 2006, ). Plus, evidence for such emission in 



([20061 . This 
, rather than 



infrared color-color ratios by Reach et al 
is likely excitation of PAHs from the ISiV 
PAHs from supernova dust. 

2. DATA USED 

2.1. Infrared Data 

Infrared data used in this paper originates from two 
Spitzer surveys: GLIMPSE using the Infrared Array 
Camera (IRAC) and MIPS GAL (survey of th e inner 

Both 
~/ < 65° 



2009) 



Galactic plane using MIPS, Carey et al 
surveys cover the Galactic coordinates lO' 
and 285° < I < 350°, |6| < 1°. 

2.1.1. MIPSGAL 

The MIPS wavelength coverage is approximately 20 to 
30 ^m at ch 1 (24 fim) and 50 to 100 /xm at ch 2 (70 
/im), thus covering a potentially rich set of diagnostic 
lines and dust emission features (see Fig. fl]). In the 24 
/im bandpass, the classic [O IV] 25. 9 ana [Fe II] 26.0 



lim tracers can be q uite strong (e.g., Oliva et al 



1999a 



Noriega- Crespo et al. 2009), and lik ewise |U 1| 63 anc 
[0 III] 88 urn (e.g., pIho et al.|[200T] ) within the 70 /xm 
bandpass. Other standard shock excitation tracers, like 
[Ne II] 12.81 or [Si II] 34.8 /im, are unfortunately not in- 
cluded within the IRAC or MIPS bandpasses. Neverthe- 
less, when dust is present, it is the dominant component 
of the IR radiation at 24 and 70 /im (see Fig. [I| . 

Although the 24 /im data delivered by the Spitzer Sci- 
ence Center in their Post-BCD products are of a very 
high quality, the MIPSGAL pipeline was designed to en- 
h ance the data a step further. The details are described 



Mizuno et al. (2008), but briefly the pipeline: (1) care- 
fully deals with saturated and non-saturated data, a key 
issue in the Galactic Plane; (2) removes a series of arti- 
facts, the most common being a column-to-column jail- 
bar striping, plus (3) that of short-duration afterimage 
latencies, (4) long-duration responsivity variations, and 
(5) background mismatches. The MIPSGAL observa- 
tions were designed to be a 24 /im survey as the high 
background levels of the Galactic plane at longer wave- 
lengths saturated the 160 /im detectors and placed the 70 
/tm at a fluence level not well characterized. Since MIPS 
takes data simultaneously at 70 and 160 /im, and that at 
70 /im at first sight looks reasonable, then a concerted 
effort has been invested by the MIPSGAL team to im- 
prove its quality. The specific steps of the 70 fiia pipeline 
are described in Paladini et al. {in preparation), but the 
main issue is to mitigate the effect of the non-linear re- 
sponse of the Ge:Ga photoconducting detectors. Their 
standard behavior leads to striping of the images, sud- 
den jumps in brightness and a large uncertainty in the 
calibration for bright sources and regions. The MIPS- 
GAL 70 /im pipeline reduces these effects and decreases 
the calibration uncertainties to a level of ~15%. This 
uncertainty on the absolute calibration of the extended 
emission is consistent with that determined b y the MIPS 
team in their calibration at 24 and 70 /im (Engelbracht 



et al. 2007 Gordon et al. 2007). The adopted uncer- 



tainty m calibration tor MIPS 24 /im is 10% . For the 
70 /zm fluxes, the main errors come from uncertainties 
in the background estimation and calibration errors (we 
adopt a conservative value of 20%). Moreover, we have 
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used the 24 /xm point source subtracted data. For details 
in the point source removal procedure see Shenoy et al. 
{in preparation). Point source sensitivities are 2 and 75 
mjy (Scr) at 24 and 70 /xm, respectively. The resolution 
of the 24 /im data is 6" while for 70 /im it is 18". 

2.1.2. GLIMPSE 



At the IRAC bands (Reach et al. 20061 in particular 
Fig. 2) one expects to iind a wealth ot atomic lines, in- 
cluding Bra (4.05 ^im), Pf/3 (4.65 ^im), [Fe II] (5.34 ^m), 
and [Ar II] (8.99 /xm), among others, plus the molecular 
S(13) through S(4) H2 rotational lines, and CO tran- 
sitio ns at 4-5 ixm. This is very much the case in IC 
443 ( [Neufeld fc Yua^[2008l [Cesarsky et "al|[l999| . We 
use primarily the 8 /im images from IRAC. Uccasion- 
ally, other channels are also used in order to constrain 
the possible emission mechanisms. The 4.5 /im chan- 
nel provi des a valua ble diagnostic of shocked molecular 
gas (e.g.. Reach et al. 2006). At these wavelengths, our 
measured tlux represents an upper limit since star con- 
tamination is present. We adopt a conservative value of 
10% for the calibration error given the vagaries in color 
correction plus measurement errors. The IRAC sampling 
is 1.2" but the 8 /zm resolution is about 2". 

2.2. Ancillary Data 

Most of the supernova remnants in Green's catalog 
(Green 2009a I are discovered using radio observations, 
ijy comparing the emission in at least two different ra- 
dio wavelengths, one can calculate the spectral index and 
thus distinguish the kind of emission produced (thermal 
or non-thermal). In the case of SNRs, we expect to find 
a spectral index which is characteristic of synchrotron 
radiation. In contrast, H II regions show a flat spec- 
trum (when optically thin) which is indicative of free- free 
emission. The SNR radio emission is associated with 
regions where shock-induced particle acceleration take 
place within the source. 

In addition to data in Green's catalog, the following 
radio data sets are used for comparison with the infrared 
images: VLA (Very Large Array) archival data at 20 and 
90 cm, available from the Mul ti- Array Galac tic Plane 
Imaging Survey (MAGPIS; Helfand et al.|2006|) websiteP 



detected (see Table[T]). Level 1 detections have a good po- 
sitional match between the infrared contours and those 
at other energies (radio or X-ray). Level 2 detections 
display infrared emission within the remnant region but 
confused with cirrus or nearby H II regions. Level 3 
detections show the presence of some infrared emission 
which does not seem to be correlated with the remnant, 
although source confusion prevents a clear assessment of 
the emission origin. Level 4 is characterized by the ab- 
sence of infrared emission within or nearby the contours 
of the radio remnant. 

Out of a sample of 121 possible radio remnants, 39 are 
level 1 detections, a detection rate of 32% compared to 
18% with GLIMPSE. A few examples of level 1 detec- 
tions are shown in Fig. |2] The MIPSGAL sensiti vity at 
24 iMo. matches very weTlthat of IRAC at 8 /tm ( Carey 
et al. 2009 Benjamin et al.||2 003). The fact that we can 
identity many more SiNRs tnan with IRAC could poten- 
tially be due to a n extinction e ffect (lower at 24 /tm than 



at 8 /im by 40%, Draine 2003). However, it is far more 



plausible that it could simply be the details of how SNRs 
emit. In the case of SNRs interacting with the ISM where 
shocks modify t he dust grain size distribution (e.g., Guir] 
lard et al.|2010 ), 24 /im continuum emission (due to Very 
Small Grains, 'V^SGs) rises while at 8 /im the emission by 
stochastically heated PAHs is not as strong. 

We calculate flux densities at 8, 24, and 70 /im for 
SNRs which show obvious infrared counterparts (level 1 
detections). The aperture size used for the flux deriva- 
tions is, in some cases, different from the one listed at 
radio wavelengths (see Table [2| . This is done in order 
to account for the different morphologies a remnant can 
have in various wavelength ranges. Furthermore, there 
are cases where the SNR is nearby other unrelated ex- 
tended Galactic sources whose contamination can lead to 
an overestimate of the infrared flux (e.g., diffuse Galac- 
tic emission and H II regions). To avoid such problems, 
we examined the MIPS images overlaid with contours 
from other wavelengths to help constrain the location 
and shape of the SNR infrared emission; however, such 
discrimination is not always well achieved. Another issue 
that arises, especially when dealing with Galactic plane 
images such as these, is the presence of unrelated infrared 
emission along the line of sight. This confusion can hin- 



VGPS (VLA Galactic Plane Survey; |Stil et al.||2006[ ]^| der th e de tection of SNRs as demonstrated by [Arendt 
at 21 cm and MOST (Mongolo Ob servatory Synthesis (jlggg^ and'Saken et al.|(|T992l) 



Telescope; Whiteoak & Green 1996[ p at about 36 cm. 

Wide band (300-10000 eVj X-ray images from the 
Chandra Supernova Remnant Catalog archive are also 
used for comparison. X-ray emission locates the shock 
heated plasma, including the important reverse shock in 
the interior region. 

3. DETECTIONS OF SNRS 

To describe the amount of (or lack of) infrared emission 
from SNRs in our sample, we use a similar c lassification 
scheme to the one adopted by Reach et al. (2006). De- 
tection levels are characterized as follows: (ly detected, 
(2) possible detection, (3) unlikely detection and (4) not 

^ http:/ /third. ucllnl.org/gps/index. html 

2 http: //www. ras.ucalgary.ca/VGPS/VGPS-data. html 

^ http:/ /www. physics. usyd.cdu.au/sifa/Main/MSC 

* http:/ /hea-www. harvard.edu/ChandraSNR/index. html 



I'lux densities ot the detected remnants were obtained 
using either a circular or elliptical aperture with a back- 
ground removed. The background was determined using 
the median value of the sky brightness around the rem- 
nant via two methods. We used an annulus around the 
source for cases where the external source contamina- 
tion in the neighborhood is negligible. However, that 
method was often not feasible due to considerable over- 
lap of sources (extended or point-like) with the periphery 
of the SNR. For those cases, the sky brightness was es- 
timated with a user-selected region characteristic of the 
local background emission. In order to quantify the mag- 
nitude of background variations, we measured the surface 
brightness on several regions around one of the faintest 
24/tm SNR, G21.8-0.6. We used aperture sizes similar to 
the one used for the remnant and found that the standard 
deviation relative to the mean background value across 
the field, for our three wavelengths are of the order of 
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12% (8 fim), 10% (24 ^m) and 20% (70 /xm). 

Another consideration for photometry measurements 
is point sourc e con tamination, especially at 8 fim (as 
discussed in § 2.1.21. Finally, since we are dealing with 
spatially extended emission, we applied corrections to the 
8 /im estimates to account for internal scattering in the 
detectors as instructed by the IRAC data handbool^ 
The 'infinite' radius aperture is most appropriate for the 
angular extent of the remnants; we multiplied the surface 
brightness measurements by 0.74. 

4. DISCUSSION 

4.1. Comparison of detections with previous infrared 

surveys 

Out of the previous 18 de tections (level 1) using 



GLIMPSE ( Reach et al 
same level with MiPSGAL 



2006), we confirm 16 at the 
i'able 1). The remaining two 



remnants are confused with the background and are re- 
ported as possible detections in our analysis. Of the other 
23 (level 1) infrared counterparts found in this work, 10 
are not in the initial sample from Green's catalog used 
by Reach et al. (2006) and the remaining 13 are mostly 
bright X-ray remnants which seem to b e brighter at 2 4 
fxm. Regarding the SNRs identified by Arendt (1989), 
all of their level 1 detections are also seen in MIPS GAL 
with the exception of G12. 0-0.1 and G54.1-I-0.3 which are 
both level 3 in the MIPS data. Comparing with |Saken| 



et al. 

(level 



1992 



, six SNRs have been detected m common 
hey also reported a level 1 detection of 8 
others, which we classify as level 2 or 3. Again, we ob- 
serve that these cases are mostly regions where there is 
a large contamination by nearby or overlapping H II re- 
gions making a clear identification challenging even at 
the superior Spitzer resolution. 

4.2. Lack of infrared signature 

Given the strong ambient emission in the inner Galac- 
tic plane, it is not surprising that many SNRs, either too 
big or too old, and having low infrared surface brightness, 
do not h ave an identifiab le infrared signature. Also, as 
noted by Williams et al. (2006b), the environs in which 
SNRs are encountered need to be sufficiently dense to 
promote collisional heating of the dust, thus making it 
distinguishable. Examples of non-detections in this sur- 
vey are G18.1-0.1 and G299.6-0.5. 



and 



4.3. Color-color Diagrams 

Table [2] shows fiux densities (at 8, 24 and 70 /im) 
color ratios for the detected remnants. Figure |3| con- 
tains a logarithmic plot of [Fg/F24] versus [F70/F24] for 
which SNRs are detected as a whole (see last column of 
Table |2|). 

Two other S NRs are plotted: for Gas A, values are 



retrieved from Hines et al. 
Noriega-Crespo et al 



more characteristic of younger remnants, such as Gas A, 
Kes 73 and Gil. 2-0.3. 

In some remnants, infrared emission is only detected in 
certain parts inside the radio/X-ray contours. Figure |4] 
shows surface brightness ratios of these localized regions. 
To assess the mechanism of the infrared emission, we also 
calculate surface brightnesses for regions of SNRs known 
to be either interacting with neighboring molecular re- 
gions (such as the BML region in 3G391) and/or that 
show molecular emission lines (such as W49B) or known 
to have ionic lines (such as 3G396). For this exercise, we 
used regions of some remnants for which the main emis- 
sion m echanism had previously been identified by [Reach] 



et al.| (2006) (see their Fig. 22). These are also identi- 
fied in the SNR figures in the appendix. Remnants with 
colors characteristic of ionic shocks seem to occupy the 
lower left part of the [I70/I24] versus [I8/I24] diagram, 
while remnants whose colors are consistent with molec- 
ular shocks and photodissociation regions are found in 
the upper right part. However, given the small number 
of examples and the fact that all three passbands (8, 24 
and 70 /im) contain dust emission and lines, we suggest 
that the [I70/I24] smd [I8/I24] color ratios are not a com- 
pletely secure method of distinguishing between different 
emission mechanisms. 

We plot some other color ratios for comparison. The 
locus for pure synchrotron is for spectral index of 0.3 to 1 
Q For the diffuse IS M, we used the color ratios obtained 
by'Gompiegne et al.'(2010) for two regions (with different 
abundances ot very small grains) in the Galactic plane at 
a longitude of approximately 59° . We have also included 
color ratios for evolved H II regions (found near Galactic 
longitudes of 35°; Paladini et al. (in preparation)). 

Are the color ratios [F70/F24] and [F8/F24] good dis- 
criminators between such very different regions? Based 
on Figure |4j the plotted colors of individual evolved H II 
regions often seem indistinguishable from SNRs using 
Spitzer colors. However, the overall trend for SNRs is 
displaced from the H II regions and co vers a broader 
range of colors. This is similar to what Arendt ( 1989 1 



found for older SNRs and compact H II regions using 
IRAS colors. 

Although line contamination is significant for radiative 
remnants, for younger ones the majority of the IR emis- 
sion should be due to dust grains. For instance, for R GW 
103 (age approximately 2000 yr; Carter et al.''1997) we 
calculate the ratio of line emission (using ISO data) ver- 
sus dust continuum within MIPS filters. We find that 
line emission contributes about 6% and 3% of the ob- 
served flux at 24 and 70 /im, respectively. 

If emission from the remnants comes primarily from 
shocked hot dust, then we can derive the dust tempera- 
ture just by assuming a simple modified blackbody emis- 
sion. If dust has a thermal spectral energy distribution 



( |2004[ ) and for IC443 from 

(in preparation). The data show 
a linear trend with a slope of 1.1 ± 0.2: remnants with 
a low [F8/F24] ratio statistically have a low [F70/F24] 
ratio. Furthermore, there seems to exist an age effect, 
with older remnants (such as IC443 and W44) showing 
higher [F70 /F24] and [Fg /F24] . The opposite seems to be 

^ http: / /ssc. spitzer. caltoch.cdu/irac/iracinstrumenthandbook/IRAC- 
Instrument-Handbook.pdf 



CX Bu{T) V 



/3 



(1) 



where /3 is the dust emissivity index (which depends on 
the dust composition, /3 ^ 1 to 2), then fiux ratios im- 
ply certain temperatures. Figure p] shows temperatures 
based on the ^3/^24] and [F24/F70] ratios. The latter 
flux ratio yields dust temperatures for our SNR sample 

^ o is the spectral index and the radio flux density is defined as 
Su oc where u is the frequency. 
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ranging between 45 and 70 K for /3 = 2 and from 50 to 
85 K for /3 = 1. This is a very crude estimate of the dust 
temperature given that we are just considering one pop- 
ulation of dust grains in temperature equihbrium when 
it is far more hkely that there are several populations 
plus non-equilibrium emission. Applying the same sim- 
ple dust model to the ^8/^24] ratio, we obtain temper- 
ature values higher than 145 K (for (3 = 2) with younger 
SNRs having color temperatures lower that older rem- 
nants. This trend is similar to what was obtained by 



Arendt] ([1989 ) using color temperatures based on the 
iKAy 12 to 25 fj,m. ratio. The infrared colors are not 
well explained by a single equilibrium dust temperature 
and atomic and molecular line emission as well as PAHs 
might significantly contribute to the infrared emission for 
SNRs. 

However, assuming that the infrared emission at 24 
and 70 ^m is cospatial, entirely due to dust, and well 
fitted by a single temperature modified blackbody, then 
SNR dust masses are given by the following equation: 



Mdust — 



Ki, B^{T dust) 



(2) 



where d is the distance to each SNRs, F^, is the flux den- 
sity and Ki, is the dust mass absorption coefficient. For 
the dust mass absorption co e fficien t, we use the diffuse 
ISM model by Li & Draine (2001) which consists of a 
mixture of amorphous silicate and carbonaceous grains. 
Distances are retrieved from Green's catalog. Results are 
reported in TablejSj Note that these derived dust masses, 
based on the 24 and 70 /im fluxes from MIPS, are overes- 
timated given the probable contamination of these fluxes 
by line emission which can be substantial in radiative 
remnants. Even in the case where all of the flux in these 
bands is due to dust, mass estimates are based on a sin- 
gle color temperature (24/70) which can be reasonable 
for younger remnants but not ideal for older remnants 
given that their IR SEP generally requires more t hat 
one dust population ( |Arendt||T989[ Saken et ari|1992[ ). 



4.4. Infrared to radio (I.4 GHz) ratio 

According to Haslam fc Osborne] (1987), the compar- 
ison between mid-infrared (bO /xm from iRAS) and ra- 
dio continuum (11 cm) emission is a good diagnostic for 
distinguishing bet ween thermal from n on-thermal radio 
emitte rs (see also (Fuerst et al. 1987) and (Broadbent 
et al.| [l989)). H II regions are shown to have infrared to 
radio ratios of the order of > 500 while rat ios for SNRs 
are thought t o be lower than 20. Moreover, Whiteoak & 
Green| ( |1996| using MOST found that SNRs have a ratio 
of infrared (60 /xm) to radio (843 MHz or 36 cm) of < 50 
while H II regions have again ratios of the order of 500 
or more. 

We us e the 1 GHz flux densities provided in Green's 
catalog (Green 2009a I. Those values were converted to 
1.4 GHz (21 crn) flux densities using the spectral indices 
quoted in the same catalog. We compare those to the 
24 and 70 /im emission from MIPSGAL. Figure [5] shows 
the spread in the ratio of the infrared to radio for the de- 
tected SNRs. Note that these ratios are often deflned by 
qmx = log(FiRA/F2icm) where Fjra and F2icm are the 
flux densities (in Jy) at specific mid-infrared wavelengths 
and in the radio. 



There seems to be a group with ratios (or high (724 
and g7o) similar to those of H II regions. In partic- 
ular, we suggest that the morphology of the remnant 
G23.6-I-0.3 (and possibly G14. 3+0.1) resembles more 
closely an H II region and so whether it is a SNR should 
be re-considered. See also its infrared colors com pare d to 
H II regions (Figs. [3]and|4|. Image in Figure 21 (and 13 1 
clearly shows a composite infrared emission, with the 24 
fim (probing hot and small particles) matching the radio, 
while the 8 and 70 /im emission do not spatially coincide 
with the 24 /im. In particular, the 8 /im appears as an 
outer layer which marks the transition between the ion- 
ized region and colder gas. The 70 fim emission is also 
seen along the same location as the 8 /xm, which is con- 
sistent with the premise that we are now probing a colder 
dust population. 

There is a strong empirical correlation between in- 
frared and r adio emission (at 1 .4 GHz) found in spi- 
ral galaxies (Helou et al. 19851. The radio emission 
in this range is thought to oe mainly due to syn- 
chrotron emission with only 10% is attribu ted to thermal 
bremsstrahlung produced by H II regions ^Mur phy|[200^ 
and references therein) . In our Galaxy and other local 
star forming galaxies, the bulk of the mid and far-infrared 
arises from warm dust associated with star forming re- 
gions. In older systems, things are more complicated 
due to cold cirrus and evolved stars (e.g., AGBs). We 
look for a similar correlation using the infrared and radio 
emission from our SNR sample. We estimate the mid- 
infrared emission as the sum of the integrated emission 
in each MIPS passband (24 and 70 fim) and obtain an 
integrated mid-inf rared flux Fmir in units of W/m^ (see 



Helou et al.||1985l ). The Fmir is an approximation (un 
derestimate) of tne mid-infrared bolometric flux (given 
that the bandpasses of 24 and 70 /xm filters do not over- 
lap). Figure shows the correlation between infrared 
and 1.4 GHz non-thermal radio flux. The slope of the 
correlation is 1.10 ± 0.13 when combining all data (plus 
IC443). Again, two distinct populations seem to exist, 
an upper main trend and a lower one. Using only the 
lower one (i.e., objects which have 70 fim to 1.4 GHz ra- 
tios closer to H II regions and high values of 924 and gro 
in Fig. [5]), we obtain a slope of 0.93 ± 0.13 while for the 
upper trend we get 0.96 ± 0.09. 

The correlation with the slope fixe d to unity leads to 
the dimensionless parameter ^mir (Helou et al. 19851 
which represents the ratio of mid-infrared to radio and 
is deflned as 



QMIR = log(^ 



F 



MIR 



3.75 X 1012 Wm-2 



- log 



F 



1.4 GHz 



Wm-2 Hz-1 



(3) 

Our results show that remnants from the lower trend 
population in Figure [6] have qmir larger than 2.06 (up to 
2.60). The other group has qum ranging from 0.06 and 
1.41. Specific q values for each remnant are displayed on 
Table Ml For reference, a previous study by BcU (2003) 
reported a median value of gxiR (i-e., the ratio of the 
total 8-1000 /im IR luminosity to the radio power) of 
2.64 for a sample of 162 galaxies. But note this would 
be higher because it is based on the total IR emission. 

Table [5] reports on the average and dispersion of the 
monochromatic qs, <Z24 and gyo parameters for the two 
trends. For comparison, a study of extragalactic VLA ra- 
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dio sources by Appleton et al. (|2004 1 obtained 524 = 0.84 
and = 2.15 whiie |ijoyle et al.] ( |2007[ ) foun d g24 = 1.39 
for so mewhat fainter galaxies. Later on, Seok et al 



( |2008[ ) using AKARI data for some LMC SNRs com- 
pared the ratio of 24 /im to radio fluxes; their correlation 
implies (724 — 0. If SNRs are the sole contributors of syn- 
chrotron emission in a star-forming galaxy, then by com- 
parison with the above (j'24's for galaxies they concluded 
that about 4 to 14% of the 24 emission in galaxies is 
due to SNRs. By doing a similar exercise using the up- 
per trend remnants, for which 1724 = 0-39 (see Table pi), 
we find that 10 to 35% of the 24 /im galactic emission 
seen would be due to the remnants. Likewise, at 70 fiui 
that contribution would be 11%. The rest would be due 
to dust heated in H II and photo-dissociation regions, as 
well as diffuse emission. 

4.5. High energy emission from SNRs 

In this section, we explore the relationship between the 
infrared and X-ray energetics of SNRs. Shocks produce 
hot plasma whose thermal energy is transferred to the 
dust gr ains via collisions and then re-emitted in the in - 



frared (Ostriker & Silk 1973 Draine 1981 Dwek 1987) 



Heating up the dust grams takes energy trom the A- 
ray gas, thus cooling it and a good tool to measure this 
transfer of energy is the ratio of infrared to X-ray power, 
usually referred to as IRX (Dwek et al. 19871, which 



compares dust cooling (gas-grain collisions] with X-ray 
cooling (continuum and lines). The gas cooling function 
increases at lower plasma temperatures given that more 
lines become available due to recombination. On the 
other hand, at high temperatures of around 10*" K, most 
of the e nergy is released throug h the bremsstrahlung con- 
tinuum ( Raymond et al.|1976 ). Assuming that line emis- 
sion in the infrared is not energetically significant com- 
pared to the X-ray lines, then most of the line cooling 
must happen in the X-ray domain. This implies that if 
the observed powers in the infrared and X-ray are com- 
parable then dust must be the essential contributor to 
the cooling in the infrared. 

Table [6] shows the values of IRX for a sample of SNRs. 
The total X-ray fluxes (with energy range from 0.3 to 10 
keV) were retrieved from the Chandra Supernova Rem- 
nant Catalog Infrared flux densities (in Jy) were mea- 
sured for regions approximately matching those in the 
X-ray. We obtain the IR flux by integrating under the 
SED using simple linear interpolation in log space plus 
i/Fv constant across the 24 /im passband. Also, note that 
the flux densities used in this analysis are not corrected 
for extinction. For our sample, IRXs range from about 
1.6 for remnant RCW103 to 240 for W44. 
These data are presented as SEDs in Figure [7j in the log- 
arithmic form XF\ vs. A which is convenient for assess- 
ing the energetic s. The radio values are obtained from 
Green's catalog (Green 2009a[ ). We also include Gas A 
and IG443 for co mparison. I'heir infrared flux densities 
are taken from Hines et ay(|2004) and Noriega-Grespo et 
al. {in preparation) while the X-ray fluxes are from the 
Chandra Supernova Remnant Gatalog and Dwek et al. 
(1987), respectively. In Figure[7|the remnants are sorted 
according to increasing age (see Table |6]). Although the 
ages for some remnants are uncertain, this Figure shows 

* http:/ /hea- www. harvard.edu/ChandraSNR/index.html 



that hot dust grain cooling (24 /im) tends to be the most 
important contributor in early phases of the remnant, 
while in the later stage the warm dust traced by the 70 
/im emission is the main contributor to plasma cooling. 
It also appears that older remnants tend to have greater 
IRXs. This can potentially be either due to an increase 
in the dust-to-gas ratio and/or a change in the size dis- 
tribution favoring smaller sizes. 

5. CONCLUSION 

We have compiled a catalog of SNRs detected within 
the MIPSGAL survey at 24 and 70 fim, with comple- 
mentary measurements at 8 /xm from the GLIMPSE sur- 
vey. In order to better assess the nature of the detected 
infrared emission, we have compared it with radio and 
X-ray data. Our main findings are the following: 

• The detection rate of SNRs given the MIPSGAL 
sensitivity is 32%, 39 out of the 121 candidates from 
Green's SNR catalog and higher than any previous 
infrared survey. 

• We find a linear trend (slope — 1.1 ± 0.2) in the 

logarithmic relationship between [F8/F24] versus 
[F70/F24]. If there is indeed an age effect, then 
the youngest SNRs will have the lowest [F8/F24] 
and [F70/F24] ratios. 

• The [I70/I24] and [I8/I24] color ratios provide a 
method of distinguishing between different emis- 
sion mechanisms. This is not completely secure 
and the color ratios of some SNRs overlap with 
those of H II regions. 

• Assuming a simple modified blackbody model (at 
24 and 70 /im), we retrieve SNRs dust tempera- 
tures which range from 45 to 70 K for a dust emis- 
sivity oi f3 — 2. 

• Using the previous color temperature (T24/70), we 
find rough estimates of dust masses ranging from 
0.02 to 2.5 M0. Note that the dust masses obtained 
here may be overestimated given the possible con- 
tribution of line emission to the MIPS fiuxes. 

• We also compare infrared fiuxes with their corre- 
sponding radio fiuxes at 1.4 GHz and find that most 
of the remnants have ratios of 70 /im to 1.4 GHz 
characteristic of SNRs, although six (about 18% 
of the detected sample) have ratios closer to those 
found for H II regions. 

• The slope of the logarithmic correlation between 
'total' mid-infrared fiux (24 and 70 /xm) and the 
1.4 GHz non-thermal radio flux is close to unity 
(1.10) as found for galaxies, qum values were cal- 
culated for each fully detected remnant and they 
range between approximately 0.06 and 2.60. 

• Whether the strong 24 /im emission is the result 
of line emission or hot dust, it is clear that there 
is a good morphological association of the 24 /im 
and X-ray features in bright X-ray remnants. The 
mechanism for the 24 /im emission for these rem- 
nants is most likely grains heated by collisions in 
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the hot plasma. The morphology of this mid- 
infrared emission is also generally distinct from the 
other infrared wavelengths which implies that the 
emission at 8 and 70 fim has a different origin. 

• We present SEDs (radio, infrared, X-ray) for a sam- 
ple of remnants and show that the energy released 
in the infrared is comparable to the cooling in the 
X-ray range. Moreover, IRX seems to increase with 
age. 
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APPENDIX 

FIRST GALACTIC QUADRANT: 10° < L < 60° 
Gl 0.5 -0.0 



First reported in Brogan et al. (2006) (hereafter, B06), this remnant was classified as a class II candidate (sources 
with the highest confidence index of detection are denominated as class I and the least as class III) . The remnant has a 
possible ASCA {Advanced Satellite for Cosmology and Astrophysics) X-ray counterpart (Sugizaki et al. 2001). Figure 
[8] shows that the 24 /im emission roughly matches the radio contours. At 70 /im, it only emits in the northern part 
and at 8 /im, the emission is confused with the background, although an arc of emission close to the X-ray source is 
well detected. This feature is also present at 24 and 70 fj,m. There is a small H II region inside the remnant which 
coincides with some of the 24 /im emission implying that the flux measurement is an upper bound. 



The infrared emission in all wavelengths is easily 
However, the brighter radio knots inside the radio 



Gil. 1+0.1 

This is a class I object from B06 with a radio shell morphology, 
confused with the background on the southern part of the remnant, 
contours show strong 24 fj,m emission (see Fig. [9]). 

Gil. 2-0. 3 

This fairly young remnant (1600 yr) contains a millisecond pulsar PSR J1811-1925 at the centre ( |Kaspi et al.|2001 ). 
It is clearly seen in th e radio as a ci rcular shell of approximately 4' diameter (spectral index of 0.6) a nd it is located at 



Rho et al. (20091 



a distance of 4.4 kpc ( Green|2009a[ and references therein). Using staring observations with Spitzer,^ ^ 

found Fe emissio n with teniperature and density representative of ejecta material in regions with high be content. 

(|2007| detected near-infrared [Fe II] and H2 filaments within and outside the radio borders of the remnant. 
While the H2 emission is most likely the result of some sort of interaction between the Shockwave and 



Koo et al. 



respective 



the progenitor wind, the [Fe II] emission is thought to be the combination of ejecta and surrounding ISM with a 
morphology characteristic of an 'asymmetric explosion'. Recently, new near-infrared spectroscopic observations of 
dense Fe knots ( |Moon et al.||2009| ) strengthened the belief that the SN explosion happened along the north-south 
direction (in Galactic coordinates). The 24um image (Fig. [lO|) reveals peaks of mid-infrared emission that coincide 



morphologically with X-ray emission (Fig. 10). Previously detected in the IRAC wavelengths, this remnant shows 8 
fiia emission in the eastern part but since it is diffuse, a clear association is not possible. The same diffuse pattern is 
seen in the 70 fj,m image. Both filaments detected in the southern rim at 8 /im have correlated emission in the 24 and 
TOjzm images, although unresolved. The IRAC colors for these filaments suggest molecular emission from shocked gas 
(Reach et al.j2006j hereafter, R06). The blob in the most eastern part of the remnant is detected both at 24 and 70 



/zm and so it is most likely due to dust. 

G12.5-h0.2 

This is a possible composite class II remnant from B06. Figure [TT] shows a striking ridge of 24 /im emission on the 
central-western part of the remnant. To the southwest it borders an IR source. Also, filamentary structure is seen 
towards the east at the boundary with a dark cloud. 

GI4.I+O.I 



This is a class II shell remnant from BOG. The 24 /tm emission has a horseshoe shape (Fig. 12 1 which roughly 



matches the radio morphology. At 8 and 70 iim only diffuse emission is seen. The remnant is surrounded by all kinds 
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TABLE 1 

Infrared detection classification levels of SNRs from 
Green's catalog 
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^ The radio sizes arc given in arc-minutes and were taken from Green's 
catalog. When two dimensions are reported they refer to the major 
and minor axis of the chipsc. 

^ IRAS survey: Arcndt 1989 , similar classification scheme as the one 
employed for this work 

IRAS survey: S aken et al. 1992. D etection levels 1 to 4 explained in 
the text. For the fSaken et al.| ( |l992| | survey, detection classification is 
as follows: Y (detected), N (not detected), P (possible detections) and 
? (not conclusive). We have translated these to the same classification 
used in this work for ready comparison (Y:l, P:2, ?:3, N:4). 

GLIMPSE survey: Reach et al. 2006. Similar classification scheme 
as the one employed for this work 
^ MIPSGAL survey: this paper 
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TABLE 2 

Characteristics of detected SNRs in the MIPSGAL 
survey'' 



1+b 


Name'=' 


SradC) 


Sphot(') 


F8(Jy) 


F24(Jy) 


F7o(Jy) 


8/24 


8/70 


24/70 


Region'^ 


10.5-0.0 


10.5 


6 


4.8 


14 lb 2 


19 it 2 


227 ± 47 


0.75 


0.06 


0.08 


whole 


11.1+0.1 


11.1 


12 X 10 


5.4 X 3.8 


12 it 2 


16 it 2 


nn'7 1 ^n 

297 ± 6U 


0.71 


0.04 


0.06 


central region 


11.2-0.3 


11.2 


4 


4.8 


4.6 it 1.1 


40 it 5 


101 it 23 


0.12 


0.05 


0.40 


whole 


12.5-1-0.2 


12.5 


6x5 


1.5 


0.9 it 0.2 


1.5 it 0.2 


7.4 it 1.7 


0.64 


0.13 


0.20 


northeastern region 


14.1-0.1 


14.1 


6x5 


1 


0.9 it 0.1 


10 it 1 


19 it 4 


0.09 


0.54 


0.05 


northwestern region 


14.3-fO.l 


14.3 


5x4 


4 


7.9 it 1.1 


30 it 4 


171 it 36 


0.26 


0.05 


0.18 


whole 


15.9-1-0.2 


15.9 


7x5 


6.5 


3.2 it 2.5 


16 it 2 


55 it 16 


0.20 


0.06 


0.30 


whole 


16.4-0.5 


16.4 


13 


6.5 X 4.2 


40 it 5 


93 it 10 


n crn i inn 

95U ± 192 


0.43 


0.04 


0.10 


northern arc 


18.6-0.2 


18.6 


6 




26 it 4 


58 it 6 


oo4 ± llo 


0.44 


0.05 


0.10 


whole 


zU.4+(J.l 


20.4 


8 


8.2 


43 it 5 


79 it 9 


1 1 rrt 1 non 

llbU ± 260 


0.55 


0.04 


0.07 


whole 


21.5-0.1 


21.5-0.1 


5 


5.2 


12 it 3 


39 it 4 


341 di 69 


0.31 


0.04 


0.11 


whole 


21.5-0.9 


21.5-0.9 


4 


1.6 


0.4 it 0.1 


0.6 it 0.1 


3.5 di 0.8 


0.63 


0.11 


0.18 


central region 


21.8-0.6 


KesD90 


20 


8.9 X 6.8 


38 it 6 


23 it 3 


A r:io 1 no 

458 ± 93 


1.6 


0.08 


0.05 


central region 


OQ G \ n o 
Zo.u-hU.o 


23.6 


10 


iO.D X O.Z 


173 it 18 


355 it 36 


5440 ± 1090 


0.49 


0.03 


0.07 


whole 


ov /linn 


Kes / o 


4 


4.8 


3.0 it 1.2 


37 it 4 


87 it 18 


0.08 


0.04 


0.42 


whole 


28.6-0.1 


28. 6A 




3.6 


7.2 it 1.0 


12 it 2 


294 it 65 


0.62 


0.03 


0.04 


nothwestern arc 




28. 6B 




2.8 


<2.0 


3.8 it 0.5 


<35 


<0.53 




0.11> 


southern filaments 




28. 6C 




4.1 


<1.9 


5.7 it 0.8 


<23 


<0.32 




0.25> 


eastern filament 


on T n o 

29.7-0.3 


Kes75S 


3 


1.4 X 0.9 


0.7 it 0.2 


4.8 it 0.5 


<29 


0.15 


0.17> 


0.03> 


southern shell 




J^.es / o W 




0.3 X 1.1 


0.6 it 0.1 


2.9 it 0.3 


5.3 it 1.1 


0.22 


0.55 


0.12 


western shell 




Kes75 






1.3 it 0.3 


7.8 it 0.8 


<34 


0.17 


0.23> 


0.04> 


whole (both sneiis) 


o 1 n 1 n n 

31.9-l-U.O 


3C 391 


7x5 


6.2 


10 it 3 


39 it 4 


395 it 81 


0.27 


0.03 


0.10 


whole 




oOoylBML 




0.8 


0.6 it 0.1 


0.9 it 0.1 


13 it 3 


0.74 


0.05 


0.06 


BML 




oOoyllN WBar 




0.7 X 0.2 


0.2 it 0.1 


0.8 it 0.1 


3.9 it 1.0 


0.25 


0.05 


0.22 


northwestern bar 


33. 6-fO.l 


Kes 79 


10 


8.4 


18 it 6 


45 it 5 


cr 1 lie 

577 ± lib 


0.40 


0.03 


0.08 


whole 




Kes79sr 




1.7 X 0.4 


<0.4 


1.2 it 0.2 


<10 


<0.36 




0.11> 


southern filaments 


34.7-0.4 


W44 


35 X 27 


38 X 28 


384 ± 44 


348 ± 39 


/^oon 1 1 onn 

6330 ± 1300 


1.10 


0.06 


0.06 


whole 


35.6-0.4 


35.6 


15 X 11 


11 


44 it 8 


36 it 5 


487 ± 99 


1.21 


0.09 


0.07 


whole 


in n n 1 


o on/^ 

3C 396 


8x6 


3.5 X 2.5 


1.7 it 0.8 


8.5 it 1.0 


34 di 11 


0.20 


0.05 


0.25 


whole 








1.4 


3.5 it 2.3 


5.26 ± 0.6 


158 di 32 


0.66 


0.02 


0.03 


southeastern plume 




o on£j Tvnx 7 




1.4 X 2.8 


0.7 it 0.2 


4.9 it 0.6 


<22 


0.14 


0.02> 


0.14> 


northwestern arc 


41.1-0.3 


o onv 
oy i 


4.5 X 2.5 


5 X 3.6 


4.3 it 0.8 


17 it 2 


34 di 14 


0.25 


0.13 


0.50 


whole 




o onvQ TP 




0.3 


< 0.1 


0.6 it 0.1 


0.6 di 0.1 


< 0.18 


< 0.17 


0.95 


ionic shock region 


43.3-0.2 


■\s r /I no 
W4yl3 


4x3 


5 


5.7 it 0.7 


77 it 8 


cr nn i i nT 

529 ± 107 


0.07 


0.01 


0.15 


whole 




W4yl3Moi 




0.6 X 0.3 


0.3 it 0.1 


0.6 it 0.1 


6.7 it 1.4 


0.52 


0.04 


0.08 


molecular shock region 




W4yl3lon 




0.6 X 0.3 


0.2 it 0.1 


3.1 it 0.3 


4.3 it 0.9 


0.09 


0.03 


0.28 


ionic shock region 


54.4-0.3 


llL'4(JJNll( 


40 


2.8 


5.9 it 0.7 


2.9 it 0.4 


38 it 8 


2.04 


0.16 


0.08 


northeastern region 




11C4U IH, 


40 


4.7 X 2.1 


3.8 it 0.5 


8.0 it 0.9 


131 it 27 


0.47 


0.03 


0.06 


north top region 


on^; Q n Q 


296 


20 X 14 


4 X 1.8 


<0.6 


1.5 it 0.2 


<11 


<0.39 




0.14> 


western arc 


on -i 1 n 1 

304.D-I-0.1 


Kes 17 


8 


8 


17 it 2 


17 it 2 


149 it 31 


1.02 


0.12 


0.11 


whole 


o 1 n o n yi 

31U.8-U.4 


KeszOAna 


12 


1.1 


0.3 it 0.1 


1.4 it 0.2 


20 it 4 


0.21 


0.02 


0.07 


northern arc 




Kes20Ase 


12 


1.2 


0.8 it 0.1 


0.8 it 0.1 


<17 


1.1 


0.05> 


0.04> 


southeastern ridge 


311.5-0.3 


311 


5 


4 


2.1 it 0.4 


3.5 it 0.5 


58 it 14 


0.60 


0.04 


0.06 


whole 


332.4-0.4 


RCW 103 


10 


10 X 8.7 


51 it 6 


101 it 11 


648 it 134 


0.50 


0.08 


0.16 


whole 




RCW103M 




1.9 X 0.4 


0.9 it 0.3 


5.9 it 0.6 


10 it 3 


0.16 


0.09 


0.58 


molecular region 


336.7-h0.5 


336 


14 X 10 


2.3 


0.6 it 0.2 


2.4 it 0.3 


6.3 it 2.5 


0.25 


0.09 


0.38 


bow 


337.2-0.7 


337 


6 


5.4 


2.4 it 0.6 


6.3 it 0.8 


13 it 5 


0.38 


0.18 


0.48 


whole 


340.6+0.3 


340 


6 


5.4 


14 it 2 


19 it 2 


193 it 40 


0.75 


0.07 


0.10 


whole 


344.7-0.1 


344 


10 


9 


4.8 it 1.9 


24 it 3 


117 it 28 


0.20 


0.04 


0.21 


whole 


345.7-0.2 


345 


6 


1.8 


<0.4 


1.7 it 0.2 


4.8 it 1.0 


<0.13 


0.05 


0.36 


southern arc 


348.5-0.0 


348 


10 


3.6 X 1.4 


1.3 it 0.2 


2.0 it 0.3 


13 it 3 


0.66 


0.10 


0.15 


central blob 


348.5-f-O.l 


CTB37A 


15 


8 


27 it 6 


27 it 4 


336 it 86 


0.89 


0.08 


0.09 


eastern shell 


348.7-f-0.3 


CTB37B 


17 


8.2 X 4 


18 it 3 


32 it 4 


401 it 82 


0.58 


0.05 


0.08 


interface 


349.7-h0.2 


349 


2.5 X 2 


2.1 


4.1 it 0.6 


32 it 4 


240 it 49 


0.13 


0.02 


0.13 


whole 



The color ratios are obtained using the flux density value for 
the partial or whole region depending on the infrared morphol- 
ogy. The radio sizes (5rad) a-i'e given in arc-minutes and were 
taken from Green's catalog. The photometry size (5phot) refers 
to the size used for obtaining the infrared fluxes. The integrated 
fluxes have not been color corrected since the mechanism for the 
infrared emission is uncertain. Also, no extinction correction has 
been applied. Furthermore, the 8 fim fluxes were multiplied by the 
extended source aperture correction factor as described in [Reach] 
[eTal. (2005a]). The tabulated errors are the combination of the 
TtT errors plus the systematic calibration errors associated to each 
wavelength. Upper flux limits are given for remnants where point 
source contamination or apparent lack of emission (3a- detection) 
is present. 

^ Name or designation used to identify remnants in Figures L3| and 
111 

Locations and areas used for photometry of subregions of rem- 
nants as identified in the figures in the appendix. 
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TABLE 3 

COLOR, TEMPERATURES AND DUST MASSES BASED OF THE 24 AND 
70 /^M FLUXES FOR SELECTED SNRS ^ 



SNRs 


^24/70 ^ 


Distance 


Mdust 


Gil. 2-0. 3 


63 


4.4 


0.034 


G15. 9+0.2 


60 


8.5 


0.081 


Kes 73 


63 


8.5 


0.11 


Kes 75 


58 


<7.5 


<0.045 


3C 391 


51 


8.5 


1.1 


Kes 79 


50 


7.8 


1.5 


W44 


48 


2.8 


2.5 


3C 396 


58 


7.7> 


0.046> 


3C 397 


65 


7.5> 


0.029> 


W49B 


54 


10 


1.6 


RCW 103 


55 


3.1 


0.18 


G337.2-0.7 


65 


2.0-9.3 


<0.018 


G349. 7+0.2 


54 


14.8 


1.6 


CasA " 


87 


3.4 


0.008 


IC443 " 


48 


0.7-2 


<0.27 



Color temperatures and dust masses are obtained using a dust 
emissivity index of /3 = 2. 

^ Approximate values for the dust temperature (K). 
Distances (kpc) talten from Green's catalog. 
Dust mass in solar masses (Mq). 

Cas A and IC443 are included for compari son. Es t imate s are 
based on the flux densities obtained by jHin es et al.| ( |2004^ and 
Noriega-Crespo et al. {in preparation), respectively. 



TABLE 4 

FmIR-RADIO RATIOS FOR SELECTED SNRS 



Name 


9MIR 


Name 


<?MIR 


GlO.5-0.0 


2.06 


G35.6-0.4 


1.37 


Gil. 2-0. 3 


0.54 


3C396 


0.06 


G14.3+0.1 


2.16 


3C397 


0.10 


G15. 9+0.2 


0.86 


W49B 


0.83 


G18.6-0.2 


2.24 


Kes 17 


0.69 


G20.4+0.1 


2.18 


G311. 5-0.3 


0.91 


G21. 5-0.1 


2.60 


RCW103 


1.07 


G23. 6+0.3 


2.43 


G337.2-0.7 


0.83 


Kes 73 


1.06 


G340.6+0.3 


1.22 


Kes 75 


0.32 


G344.7-0.1 


1.41 


3C391 


0.87 


G349.7+0.2 


0.76 


Kes 79 


1.05 


CasA 


-1.05 


W44 


1.02 


IC443 ^ 


0.53 



^ Va lues for infrared an d radio emission from Cas A were taken 
from |Hines et al"] | |2004| | . 

^ Values for infrared emission from IC 443 were taken from Noriega- 
Crespo et al. [in preparation). 



TABLE 5 

QlRx VALUES FOR DIFFERENT WAVELENGTHS 



-JIRa 


All sample 


Upper trend 


Lower trend 


'?24 


0.71 ± 0.42 


0.39 + 0.11 


1.68 + 0.05 


970 


1.55 + 0.60 


1.17 + 0.19 


2.70 + 0.04 


Qmir 


1.19 + 0.53 


0.83 + 0.15 


2.28 + 0.04 



524 and 570 parameters are monochromatic. They are defined 
as iJiR;^ = log(FiR^/F2icm) where Fir^ and F21 cm are the flux 
densities (in Jy) at specific mid-infrared wavelengths and in the 
radio. Qmir. is a dimensionless parameter which represents the 
ratio of mid-infrared to radio (it is defined in Eq. Isl. Individual 
Qmir values for each remnant are presented in the Table |4] 
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TABLE 6 

Infrared to X-ray ratios for selected SNRs 



JN anie 


T'^+cil IV/TTR T?liiv 

(10 ^ erg cm s ^) 


A-ray 
(10- 


i^l"X(o 3_io KgV) 
-Q —2 —1 \ 

" erg cm s ) 




A (Tilt) 

Age 


m 1 2-n S 


11 




4.0 


2.8 


1 finn f 1 1 


G15. 9+0.2 


5.3 




0.96 


5.6 


500-2400 (2) 


Kes73 


9.9 




1.8 


5.6 


500-1000 (3) 


Kes75 


2.7 




0.21 


13 


884 < /> 2400 (4) 


3C391 


20 




0.62 


32 


4000 (5) 


Kes79 


16 




0.17 


93 


6000 (6) 


W44^ 


2.0 




0.093 


240 


20000 (7) 


3C396 


3.0 




0.19 


16 


7100 (8) 


3C397 


8.1 




1.3 


6.2 


5300 (9) 


W49B 


30 




6.0 


5.1 


1000-4000 (10) 


RCW103 


27 




17 


1.6 


2000 (11) 


G337.2-0.7 


1.6 




0.93 


1.7 


750-3500 (12) 


G349. 7+0.2 


14 




0.33 


42 


2800 (13) 



X-ray fluxes were retrieved from the Chandra Supernova Rem- 
nant Catalog. 

^ Age is given in ye ars and the values were take n froi n the fol- 
lowing literature: (1) [Clark fc Stephcnsonl ( |1977^ ; (2) [Reynolds] 
[eTal. (2006); (3) |TianT^°F eahy (2UUSa) ; (4 ) upper liniit lor th e 
spm-down age of the associated pulsar Livingstone et al.' f2006|/ 
based on ionization time-scales Morton et al. (2UU7I; (5) (J hen et al 
1^04); (6) Sun et al. (2004) ; (7) Wolszczan ct al. (1991); (8:tlaF 
ms & Slanc (1999); (9) Safi-Harb ct al. (2005); (10) e.g. Hwa^ 
e?al. (.2000 J: (11), Carter et al. 119971; ( 12) .Kakowski et al. (,200W 
^ ' ' 



Infrared and X-ray flux values are for a central rectangular region 
of W44. 




J I L 



20 



30 40 50 

Wavelength(/2m) 



60 70 



50 90 



Fig. 1. — A sketch of a dusty SNR spectrum made with a blackbody at 60 K and superimposed the most important fine structure and 
ionic emission lines (mimicking ISO Short and Long wavelength spectrometer observations of SNR R,CW103). Emission lines within the 
wavelength range of the MIPS filters are annotated. Overplotted (dashed lines) are the normalized response functions of the MIPS filters 
at 24 and 70 fim. This shows that the contribution of line emission to the MIPS filters is relatively small when a dust continuum is present. 
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Fig. 2.— Examples of level 1 detections at 24 /^m. Prom left to right, top to bottom: Gll.2-0.3, G27.4+0.0 (Kes 73), G31.9+0.0 (30 
391), G43.3-0.2 (W49B), G332.4-0.4 (ROW 103) and G349.7-0.2. The greyscale is Unear with the ranges in MJy/sr as displayed in the 
intensity bars below e3x:h. image. The angular extents of the images differ. 
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Fig. 3. — Color-color plot of the SNRs in our sample detected and measured as a whole. For comparison, we have also plotted the location 
of Cas A and IC443 in this diagram. We find that this sample is well fitted by a slope of f.f it 0.2 in this logarithmic representation, 
showing that 8 and 70 [iva emission are related roughly by a constant while they change dramatically with respect to 24 /xm. Vertical lines 
show the range of expected temperatures for a modified blackbody using only the flux measured at 24 and 70 /xm while for the horizontal 
lines, the temperatures are calculated using the flux at 8 and 24 ^lra. Dotted and dashed lines correspond to a /3 of f or 2, respectively. If 
the IR emission at longer wavelengths is mainly due to dust then the range of temperatures for our sample of SNRs is roughly between 50 
and 85 K for /3=1 and between 45 and 70 K for ^=2. 



of infrared filaments which without the radio contours would be hard to distinguish from the overall structure. 

G113+0.1 

This is another class II remnant from B06. While the 8 /im IRAC image shows emission localized in the eastern 
part, the 70 /im emission fills the radio contours uniforml y. T he 24 /im emission also fills the contours with two peaks 
located north and south of the remnant's centre (see Fig. 13 1. 



G15. 9+0.2 



Ch andra observations showed that this remnant is young (< 2400 yr) with an X-ray point source (Reynolds et al. 
20061. The obtained X-ray spectrum implied high element abundances characteristic of ejected material. I'he infrared 
detections are presented in Figure [Til This remnant was not detected in either the GLIMPSE or IRAS surveys. The 
24 /im emission is bright in the southern-most part. There are also traces of emission in the north rim. At 8 /im, 
the emission is very faint and mainly seen around the south-eastern radio contours. The 70 /im image shows diffuse 
emission extending along the south-eastern part; nevertheless, it is unclear if this is associated with the remnant. Also, 
in the sa me area, the remn ant seems to be encountering interstellar material, although no maser emission has been 
reported ( [Green et ar||1997[ ). 



G16.4-0.5 

B06 classified this as a class II partial shell. The 24 /im emission has a shape similar to a fork with the strongest 
emission coming from its prongs. The western arc feature does not seem to have a st rong counterpart at the other 
infrared wavelengths, unlike the rest of the structure which shows up diffusely (Fig. 15). 
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Fig. 4. — Color-color plot based on surface brightnesses of localized regions in SNRs in our sample. Ratios for SNRs known to have 
molecular interactions, atomic fine-structure line emission or PDRs colors are also shown. ICE and ICW refer to IC443, eastern (strong 
[O I] emission at 63 /im) and western (interacting with a nearby H II region) regions, respectively. Colors of pure infrared synchrotron 
emission are plotted as a straight line. Color ratios for the diffuse ISM |Compiegnc et al. 2010) and evolved H II regions (Paladini et al. (in 
preparation)) are also included. Remnants with upper flux limits are represented m italic, b'or comparison, we plot the slope (dotted-dashed 
line) of 1.1 it 0.2 found for the sample of detected remnants and measured as a whole (see Fig. [sj. Also depicted is a dashed line showing 
a different trend for evolved HII regions. 



G 18. 6- 0.2 

Another partial shell class II remnant from BG06 is shown in Figure [16] There is a good correspondence between the 
radio contours and the infrared emission, in particular, for 24 /im. At 70 /xm the strongest infrared emission matches 
the brightest radio peaks and, at 8 /xm, there is diffuse emission which seems to occupy those regions as well, but no 
clear association can be established based on this filter alone. 



G20.i+0.1 

Figure [17] shows a shell class I remnant from B06. The 24 and 70 /.tm emission fill the radio contours. The 
northwestern radio shell is well matched by the 24 and 70 /im emission and it is surrounded by 8 /im filaments. 



G21.5-0.9 

Tian & Leahy ( 2008b') suggested that this remn ant is at 4.8 Kpc ba sed on HI and CO observations and it has an 



associated pulsar (PSR J1833-1034) at its centre (Camilo et al. 20061. Figure 18 shows the infrared emission at 24 
and 70 ^m. associated with the central X-ray contours. 'I'here is no distinctive association between the 8 /im diffuse 
emission and the Chandra image. 



G 2 1.5- 0.1 

This is a partial shell remnant identified in BG06, with a possible ASCA X-ray counterpart (Sugizaki et al. 112001). 
Centred at coordinates 21.6 and —0.1, a ring with 0.5' radius accounts for most of the 24 /im emission (Fig. 
southern filament which coincides with a bright radio ridge is seen in all of the infrared wavelengths. Below tF 
is a filament of 8 //m immediately outside the radio contours. 
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Fig. 5. — Ratio of MIPS 70 /im and radio continuum at 1.4 GHz (570 = log(F7o/im/i^i.4GHz)) versus ratio between 24 fira and radio 
again ((J24)- There is a considerable range in the positions of the SNRs. The ratios are correlated and are well fitted by a slope of 1.2. 
There seems to exist a group of catalogued sources which have ratios more characteristic of H II regions (upper population in the plot). 
The infrared counterpart s in the lower part of th e plot show ratios of infrared (70 /im) to radio comparable to what was found for SNRs 
in previous studies (e.g., [Broadbent et ah] | |1989[ l) which suggested ratios of infrared (60 /.tm) to 2.7 GHz lower than 20. The location of 
. Gas A tails ott the plot and below the trend with ratios (0.1, 0.05). 



IG443 is also shown. 



G 21. 8-0. 6 or Kes 69 

At a distance of 5.5 Kpc ( Tian fc Leahy|[2"008b ) , Kes 69 was detected in the GLIMPSE survey and shocked hnes are 
most Hkely responsible for the majority of the emission seen in the south part of the radio contours (R06) . The 8 and 
24 /j,m images (Fig. 20) show filaments that might be associated with the remnant centre {Kes69C), though there are 
similar faint structures to the south beyond the radio contours. There is prominent mid-infrared emission north of the 
remnant due to a bright H II region. 

G23.6+0.3 



This remnant has a radio spectral index of -0.34 (Shaver & Goss 19701. Its odd elongated radio shape is seen to 
spatially coincide with the 24 /im emission (Fig. 21 H I'he 8 and 70 /zm emissions are more prominently displaced 
towards the south and west of the radio boundaries. No maser association has been reported. 

G27. 4+0.0 or Kes 13 

This remnant has a 4' diameter shell-like structure in the radio 
be as far as 9.8 kpc and be 500-1000 yr old 
the central compact source (IE 1841-045) 
seen at 24 /^m (Fig 



Tian & Leahy ( 2008a ) reported that this SNR can 



Vasisht & Gotthelf (1997) discovered a low period (12s) X-ray signal from 



the northwestern ridge 



he yiNK is not detected at the IRAC wavelengths but it is spectacularly 
22|). This mid-infrared emission traces the X-ray contours extremely well, with both peaking in 



G 28. 6- 0.1 



This region is highly contaminated by thermal emission and therefore extr emely bright in infrared wavelengths. The 
three areas chosen to calculate infrared fluxes are identified in Figure [23] by A, B and C. These sources showed a 



16 



—3 



N 
X 

o 



2 



O 







-1 



n 1 1 1 1 1 1 1 1 1 1 r 



n 1 1 r 



CasA 



W4.4 ■ ' 



IC443 



3C396 
Kes75 

^ ^1 1 



3C397 



VV49e 
^riof ''RCW103-' 
11.2 349 ^^^Mes79.--- 



Kesl 7 



35.6 



.-25.6 



J 5:9 



Kes73 



344 



.340 



20.4 



18.6- 



10.5 
14.3 .. 



21 .5-0.1 



_] I £l L. 



_] I I I I I I I L. 



_l I I L. 



_] I I L. 



12.5 



2.0 -11.5 

Log f^M\R) (W m 



1 1 .0 

2 



10.5 



0.0 



Fig. 6. — Correlation between radio and infrared emission. For the entire population, the slope of the correlation is 1.10 (dotted line). 
There seem to exist two trends. The upper population is more characteristic of SNRs while the lower population looks more like H II 
regions. The latter are the same objects with the high values of 524 and 570 in Figure |5] The slopes for the upper and lower populations 
are 0.96 (dashed line) and 0.93 (dot-dashed line), respectively, not significantly different. Again, IC443 and Cas A are also included for 
comparison. Cas A stands out due to its strong radio synchrotron emission. 



non-thermal spectra in the radio ( Helfand et aL|1989 ) and 20 cm flux densities of 0.73, 1.05 and 1.05 Jy, respectively. 
The radio morphology of all sources trace closely the 24 /zm shape. Sources B and C do not show any 8 and 70 /Ltm 
emission and are only well detected at 24 fim. 

G 29. 7- 0.3 or Kes 75 

This composite remnant has a n associated pu lsar (PSR J1846-0258; Gotthelf et al.|[2C)00 l) and is thought to have an 
upper distance limit of 7.5 kpc (Leahy &: Tian| [2008 ). This remnant was not detected in GLIMPSE or IRAS surveys 
but was well identified using targeted MIPS OTOervations by Morton et al. (2007) where most of the emission at 24 /im 
originates from two bright shells located south and west of the remnant's centre. Figure [24] shows a clear detection at 
24 /xm which traces accurately the brighter X-ray regions with the central X-ray peak being the only exception; there 
is no 8 /zm emission there and the 70 /.tm emission appears to be non-existent too. There is also no significant infrared 
emission at 8 or 70 /im associated with the bright 24 /im features. Unrelated 70 /im emission shows up diffusely in the 
eastern part of the remnant which overestimates the 70/24 color ratio measured for this region (Kes75S). 

G3 1.9+0.0 or 3C 391 

Green's SNR catalog describes this object as a shell-type remnant with a sp ectral index of a bout 0.49. Previous 



infrared spectral observations with ISO revealed the existence of ionic lines (Reach & Rho 2000), in particular [Fe II] 
at 26 /im and [S III at 35 /im towards the direction of the BML (Broad Molecular Line j region at Galactic coordinates 
31.84, 0.03 (Fig. 
in the radio pea' 
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Besides the previous elements, [O IV] at 26 /im was also detected but found to be 3 times brighter 
3ar at Galactic coordinates 31.8 7, 0.06) than in th e BML region. Right above this radio peak, the 
remnant is encountering a nearby molecular cloud ( Wilner et al.|1998| ^ All of the abo ve line emitters cont aminate the 
MIPS channels and therefore contribute to fiuxes reported in I'abie [2| Furthermore, Reach et al. (2002) also found 
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Fig. 7. — SEDs of Chandra SNRs ordered according to approximate remnant's age (increasing left to right, top to bottom). For the 
radio, the synchrotron spectrum is shown by the dotted line for the appropriate spectral index, anchored in the observed range at 1 GHz 
(triangles). The circles indicate the power in the infrared at 8, 24 and 70 /im. In the X-ray, the band (0.3-10 keV) of the Chandra images 
is used (squares). Horizontal lines represent the width of the bandpass. We also include Cas A and IC443 for comparison. Their infrared 
flux densities are taken fromlHines et al. (2004) and Noriega-Crespo et al. [in preparation) while the X-ray fluxes are from the Chandra 
Supernova Remnant Catalog and Uwck ct al. (1987), respectively. Again, note that for W44 the infrared and X-ray values apply only to the 
central region while the radio estimate is tor the whole remnant. The integrals over the infrared and X-ray range reveal their contributions 
to cooling the shocked plasma. 

clumps of H2 emission in the BML which was again confir med with IRAC color ratios by R06. Two 1720 MHz OH 
masers were detected within the remnant ( Frail et al.||1996 ). IRAC images show emission associated with the western 
maser where the previously mentioned BML is located. Mid-infrared 24 and 70 fim emission from MIPSGAL is clearly 
seen in this region as well (Fig. 25 ) . Previous IRAC observations revealed a rim of [Fe II] emission which matches the 



brightest radio emission in the north part of the remnant. The same feature is seen in the 24 /zm image which suggests 
that the bulk of emission at this location can be due to [Fe II] at 26 /zm. R06 also found that the infrared emission 
in the middle knots (east and west at the 'waist'; the one in the west has the same location as the BML) is produced 
by shocked molecular gas. These knots are well detected at 24 /xm (with some likely contribution from H2 S(0) at 28 
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/im) and are especially strong at 70 /ini. The southeastern rim structure (at 24 fim) coincides with an enhancement 
in the radio and the second maser. 

Remnant 3C391 is characterized by a faint outer rim that can be seen both in the radio at 20 cm and in the infrared 
at 24 iim. This is the collisionless shock where the blast wave encounters the ISM. Although the radio structure is 
spatially correlated with the faint mid-infrared, a spectral index analysis shows that this 24 /im emission is too strong 
to have only a synchrotron origin. 



G33. 6+0.1 or Kes 79 



This a 'mixed- morphology' remnant with a centrally compact X-ray source (Seward et al. 2003). It has an incomplete 
radio shell of 10' in size and it is X-ray centre filled. The southern filaments seen in the Chandra image are well detected 
and have a similar morphology in the 24 /xm image (Fig. 26). Those are located in the same regi on as a dark cloud, 
where previous CO and IICO+ observations suggested an interplay with a nearby molecular cloud d Green fc Dewdney| 
1992). These authors also noted that t his must be an evolved remnant given its large linear size, furthermore, using 



[T9S 
Ilhi 



iiandra observations [Sun et al. ( 2004[ ) did not encounter evidence of ejecta material, reinforcing the conclusion that 
this must be an older remnant (around 6000 yr old) . The 24 fim emission fills and traces the brightest X-ray contours 
(so called inner-shell) in the central part of the remnant while 8 /zm diffuse emission is present in the northern part of 
the remnant, but not likely to be associated. Likewise, at 70 /im it appears that the main contribution to the flux is 
probably the IRAS source 18501-1-0038 (represented by a white circle in Fig. 26). 



G34.7-0.4 or W44 

This evolved object is classified as a 'mixed-morphology' remnant (Rho & Petre 1998) due to its radio shell with 
inner X-ray filled emission. It is mainly detected at 70 fim possibly indicating that it has had time to evolve (see 
Fig. 27 ). This is in agreement with the inferred age of about 20000 yr from the characteristic age of the associated 
pulsar (Wolszczan et al. 1991). In the 70 fim image, the western filaments trace well the radio contours. In all the 
infrared wavelengths, there is emission associated to the northern region. This is one of the most stunning detections 
in GLIMPSE data where R06 found tha t emission at 4. 5 /j,m (which is mostly from shocked molecular gas) matches 



well the radio contours. Previously, Reach et al. (2005b I also detected near-infrared shocked II2 and millimeter wave 
molecular lines, which revealed the interaction of the re ninant with a CMC. Moreover, the existence of extended maser 
emission has also b een reported (Hewitt et al. 20081 inside of the remnant. Recent Spitzer spectral line analysis 
( [Neufeld et al]|2007[ ) revealed the emission of [be llj around 26 /zm and some H2 at 28 /zm which can be contributors 
to the taint emission seen in the MIPSGAL 24 /im image. At 8 /im the object is easily confused with the diffuse 
background. 



G35.6-0.A 



This object was recently re-classified as a SNR ( Green|[2009b ). It has a spectr 

F 15' X 11'. IRAS emission has not been previously reported. Its emission at 8 and 70 /im can 



tral index of 0.47. The radio emission 

is peaked, with size of 

be more easily confused with the background material and there seems to be no convincing association with the radio 
contours (Fig. 28). At all wavelengths, there is emission associated with a gamma-ray source HESS 1858-1-020 and a 



narrow ridge (jets?) can be seen at 24 /im. 



G39.2-0.3 or 3G396 



Several filaments are observed at 24 /im (Fig. 29). Specific ally, the to p interior part of the remnant is only observed 
in this wavelength. Also, recent near-infrared observations (Lee et al. 2009| ) revealed the presence of [Fe II] and H2 
in the northwestern part of the shell matching the radio emission, i'he 24 /im image shows the exact same features 
which suggests that some emission is probably due to lines in the MIPSGAL band. Other filaments near the centrally 
projected hot point source (IRAS 19(]17-}-0522) are seen at 8, 24 and 70 /xm. Previously, R06 had argued that those 
filaments had colors consistent with PDRs . In that case, the suit e of Spitzer pas sbands is probing the dust emission 
from PAHs to big grains. The tail of radio emission mentioned in Patnaik et al.| (199 0) is partially traced by 24 and 
70 /im with some spatial mismatch. There is a bright filament in the south part of this tail, which is detected in all 
of the infrared wavelengths. Moreover, in the eastern part of the tail, a filament of infrared emission also seems to 
match exactly the radio contour. The infrared emission that is in between the tail and the remnant is most likely to 
be caused by thermal dust grains. 



G4I.I-O.3 or 3C 397 

This is another shell- type morphology remnant with spectral index of 0.48 ( Green'2009a') . Its uncommon rectangular 



shape on the upper edge suggests in teraction with the ISM, specifically a dense molecular cloud inferred through CO 



observations (Safi-Harb et al. 2005|) 
5300 yr old and is now starting t 
part of the remnant (Fig. 
X-ray counterpart. The G 



emission from ionic lines (R06). 



Based on Chandra data, the same authors propose that the remnant is about 
le radiative stage. The 8 and 70 /im images show some emission in the north 
30). The 24 /im image reveals rich filamentary structure that nicely resembles the entire 
IMPSE survey reported a faint detection in regions where IRAC colors suggested shocked 
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G 43. 3- 0.2 or W49B 

W49B is characterized by an odd 'barrel'-shaped morphology in the radio (Moffett &: Re ynolds|1994 ). R06 reported 



that most of the emission at the IRAC wavelengths is consistent with shocked ionic and molecular emission lines. 
Later on, near-infrared observations revealed that the [Fe II] morphology, at 1.64 /zm, is similar to the radio structure 
and the 2.12 nm H2 can be found p rimarily in the south (right below the [Fe II] ridge) and in the northwestern part of 
the remnant ( Keohane et al. 2007 ) . With the help of Chandra data, these authors also reported an overabundance of 
heavy elements"thu s suggesting that the X-ray emission is dominated by the ejected material as found previously by 
Hwang et al. (20001 using ASCA data. Both MIPS channels (see Fig. 31) show interesting filamentary structure and 
are brighter in the same locations, thus indicating that there is a strong dust continuum contribution as the primary 
emission mechanism in the far-infrared. However, note that the 24 fim image seems to trace well the ionic shocked 
component within the X-ray contours. 



G54.4-O.3 

GLIMPSE data (R06) revealed two regions with the same color ratios as photodissociation regions whose infrared 
filaments spatially match the radio profile. Those are also seen as filaments in the 24 /im and diffuse emission at 70 
fim (see Fig. 32 1 . Fluxes for those two regions are reported in Table |2] 



FOURTH GALACTIC QUADRANT: 295°<L<350° 
G296. 8-0.3 



This rem nant with an elongated radio shell is at 9.6 kpc and has an age between (2 — 10) x 10 yr (Gaensler 
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et al.|[T998 ). These authors also concluded that its radio morphology might have been shaped either due to mass loss 
from winds or interaction with the ISM. There is an arc of bright 24 /im emission which matches the strongest r adio 
synchrotron emission in the northwest region. No IRAS or GLIMPSE detection has been reported. As seen in Fig. 
there is some diffuse 70 /im emission around it that is unlikely to be associated with the remnant. 

G 304. 6+ 0.1 or Kes 17 

Infrared emission from the northwestern part was detected in the GLIMPSE data and attributed to shocked molecular 
gas (R06). Those filaments are also present in the 24 fjjn image (Fig. 34 1. Moreover, at that same wavelength, there is 



emission along the south which spatially coincides with the radio ridge. At 70 /im, there is just filled emission around 
the northwestern shell. 

G3 10.8-0.4 or Kes 20 A 



MO ST observations show a well-defined eastern radio shell with a weak western counterpart (Whiteoak & Green 
19961. Diffuse thermal emission in the western counterpart of the remnant has led to previous confusion regarding 
the thermal or non-thermal origin of this SNR. Following the radio contours to the south west of the bright compact 
eastern radio peak is an arc of 24 fMn emission (Fig 



distinguishable at 8 /im. 
shell. 



35|. This feature is also present at 70 /im but it is not clearly 
There are also filaments of diSuse infrared emission that follow the north-south eastern radio 



G3 11.5-0.3 



Whiteoak & Green ( 1996 1 first identified this as a shell type SNR. With GLIMPSE data, R06 detected an incomplete 
shell (brightest in the western part). Color-color ratios led to the conclusion that infrared emission was due to shocked 



molecular gas. Both MIPS images show the same approximate morphology following a circular shell (Fig. 36 ) 



G332.4-0.4 or RCW 103 



This young remnant (age about 2000 yr) has a shell-like radio morphology and is at a distance of 3.1 kpc (|Gree^ 
2009a and references herein). It has a soft X-ray point source near the centre of the remnant (Tuohy &: Garmire 
TUSOj). The 24 /xm image (Fig. 37) shows extensive diffuse and filamentary emission related to the remnant, spatially 
correlated with the radio. Numerous small filaments trace the outer borders of the synchrotron emission. There is 
a dark arc in the centre left of the remnant where the IR and X-ray shells seem to be incomplete. For the southern 
bright rim, R06 found colors indicating shocked mo l ecular gas along with some ionic emission, which corroborated the 



reported detection of ionized species by jOliva et al. (1999b). Around 24 /im, the ISO spectrum showed strong emission 
mainly from [Fe II] as well as [O IV]. These authors also reported that heavy elements are in the gas phase given their 
solar-like abundances and therefore most of the infrared emission in that rim must be the result of fine-structure lines. 
The 24 /im morphology is very similar to the X-ray emission, as seen in other young remnants. 

G336.7-t-0.5 



There is an arc in the 24 /im image which matches roughly a radio counterpart emission peak (see Fig. 38 ) . At 
8 and 70 /tm, there is some diffuse emission in the south which coincides with a brighter radio region but it seems 
unlikely that they are associated. The flattened indented outer radio contours to the south west coincide with with a 
dark cloud, possibly the locus of interaction of the SNR with the ISM. 
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G337. 2-0.7 



The rad io contours from MOST observations indicate that this is a compact circular shell remnant (Whiteoak & 
[Green 1996). It has not been detected in the IRAC bands or IRAS surveys. The remnant is more complex m the J^-ray 
and mid-infrared (Fig. 39 ) , with the 24 /im emission composed of several filaments forming two distorted apparent 
shells. The interior one is almost complete and is located in the western part. The other is visible as two long 
filaments mainly in the north-south direction in the eastern part of the remnant. Both relate closely to the contours 
of the Chandra X-ray emission. There is some 70 /xm diffuse emission in the north-eastern part but unlikely to be 
related with the remnant. 

G 3 40. 6 +0.3 



This is a circular shell remnant in the radio located at 15 kpc (Kothes &; Dougherty 2007). Filaments of 24 /xm 
emission follow the radio contours especially in the south west (see I'lg. |40[). Also, emission from two projected IRAS 
point sources, 16441-4427 and 16441-4429, situated in the north region ofthe SNR at opposing sides, enhance the flux 
considerably at the MIPSGAL wavelengths; note that in the Spitzer bands there is considerable extended emission at 
these peaks. 

G3U.7-0.1 



Dubner et al. ( 1993 1 reported an angular size of 33 pc in the radio which agrees wi th the e s timat e of a middle-aged 
remant as based on the X-ray plasma temperature found using ASCA data Yamauchi et al. (2005). Figure 41 shows 
that the emission at all three infrared wavelengths is brightest at the central radio peak. Previous GLIMPSF analysis 
showed that this structure had colors consistent with shocked ionized gas (R06). Although now clearly diffuse emission, 
this peak was catalogued as IRAS 17003-4136. The 24 /im image also shows a few striking filaments associated with 
the outer radio boundaries, especially where compressed at the northwest. 

G345. 7-0.2 



A rim of emission at 24 and 70 /im is detected in the south-eastern part of the rem nant (Fig. 42 1 , poss ibly shocked 



dust from interaction with the adjacent dark cloud. No maser association is reported (Green et al. 



1997). 



G348. 5-0.0 

This SNR is independent from G348. 5+0.1 (see below), although its radio morphology looks like a jet of material 
emanating from CTB 37A ( Kassini et al.|[T99T ) . GLIMPSE data revealed that the arc of associated infrared emission 
is mainly dominated by ionic shocks (R06j. l<'igure 43 shows the same structure at 24 and 70 /im, the 24 fj,m emission 
being the strongest of all. 

G348. 5+0.1 or CTB 37 A 

GLIMPSE analysis supp orts the idea that this remnant is interacting with dense interstellar gas, first indicated 
by the presence of masers (Frail et al. 1996) within the remnant. Reynoso fc Mangum ( 2000^ found CO clouds at 
approximately the same velocities and locations of the OH 1720 MHz masers. li igure nSl sEows that the MIPSGAL 
images have the same diffuse emission seen at the IRAC wavelengths including the eastern patch of infrared emission 
which is thought to be physically associated with a maser (and whose IRAC colors suggest emission from shocked 
molecular gas as noted by R06). At 24 /im there are also southern filaments which closely trace the radio boundaries. 
Bright infrared emission on the northwest fills the faint radio emission part of the remnant and its IRAC colors are 
similar to photodissociation regions (R06). 

G348.7+0.3 or CTB 37B 

In this very confused field. Figure [44| shows that the radio peak from MOST contours has an associated mid- infrared 
peak at 24 /im. In addition there are some filaments within the remnant that seem to be the result of interaction with 
the ISM. Chandra follow-up obse rvations of a Gamma-ray source (HESS J1713-381) thought to be associated with 
CTB 37B (Aharonian et al. 2008 and references therein) also detected X-ray thermal emission close to the remnant 
centre whose temperature suggested a middle-aged SNR (about 5000 yr). 

G 34 9. 7 +0.2 



OH maser emission was found in the r emnant's interior by F rail et al. (1996) and evidence of interaction with a 
molecular cloud was later strengthened by Reynoso & Mangum ( |200l[ ) with CO observations. With extre mely s imilar 



2005 ). The 
I'he infrared is strongly 



X-ray and radio morphology, its corresponding brightness in both regimes is also high (Lazendic et al 
infrared emission is present with a similar morphology at all three wavelengths (Fig. 45 F 

peaked within the X-ray contours, but not closely related, except for the 24 /im emission^'revious GLIMPSE analysis 
showed infrared emission in all of the IRAC bands (R06) 



INDIVIDUAL IMAGES 



Images of the infrared emission for each detected remnant follow. The order is the same for all the figures and 
consists of the black and white 24 /im image (colorbar in units of MJy/sr) on the left and the 3-color image on the 
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right, with 8, 24 and 70 /Ltm being blue, green and r ed, respectively. Overlaid in each image are contours either from 
radio or X-ray observations when available (see { 2.2 1. Images are all in Galactic coordinates, and our remarks on north 
(up), east (left), etc., are with respect to this system. Also indicated is the angular scale which can be compared with 
the few arcminute resolution of IRAS. 




36 40 



48 62 



Fig. 8. — SNR GIO. 5-0.0. Contours from VLA 20 cm observations: levels are 13, 17, 21 and 24 mjy/beam. The cross represents the 
location of the associated X-ray source (AX J180902-1948). 
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Fig. 9.— SNR Gll.H-0.1. Contours from VLA 20 cm observations: levels are 14, 18, 22 and 26 mJy/beam. 
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Fig. 11. — SNR G12.5+0.2. Positions of a maser (diamond) and nearby infrared IRAS (circle) source are overlaid. Also shown is a circle 
showing the approximate size and location of the remnant in the radio as seen in Fig. 2 from B06. 



23 







IZ ..... . ^^^^M^^^M^^^M^^ 

36 37 39 44 50 57 67 7B 91 105 121 

Fig. 13.— SNR G14.3+0.1. Contours from VLA 20 cm observations; levels are 16, 20, 24 and 28 mjy/beam. 
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Fig. 14.— SNR G15.9+0.2. Contours from VLA 20 cm observations: 
the location of an X-ray source (CXOU J181852. 0-150213) 
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Fig. 18. — SNR G21.5-0.9. Contours from Chandra observations: levels are 2, 3, 6.5, 12 and 20 X 10 ^ photons/cm^/sec/pixel. The 
cross locates PSR J1833-1034. 
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Fig. 19. — SNR, G21. 5-0.1. Contours from VLA 20 cm observations: levels are from 7 to 10, in steps of 1 mjy/beam. 
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Fig. 21. — SNR G23.6+0.3. Contours from VLA 20 cm observations; levels are from 10 to 22, in steps of 3 mjy/beam. 
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Fig. 24. — SNR G29.7-0.3. Contours from Chandra observations: levels are 1.3, 1.8, 3.2, 5.5 and 8.7 X 10"'^ photons/cm^ /sec/pixel. The 
cross marks the location of the pulsar AX J1846. 4-0258. Regions used for partial photometry are also indicated in the figure (in black). 




Fig. 26. — SNR G33. 6+0.1. Contours from Chandra observations: levels are 1.5, 2, 3.6, 6.3 and 10 X 10~^ photons/cm^/sec/pixel. The 
cross marks the location of PSR J 1852+0040 and the circle represents an IRAS point source. 
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Fig. 28. — SNR G35.6-0.4. Contours from VGPS 21 cm observations: levels are from 18 to 30, in steps of 3 Jy/beam. The cross marks 
the location of 7-ray source J 1858+020. 
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Fig. 29. — SNR G39.2-0.3. Contours from 20cm VLA observations: levels are 1, 4, 12, 26 and 45 mjy/beam. The cross represents the 
location of the X-ray source J190406. 5+052646. Regions used for partial photometry are also indicated in the figure (in white). 
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Fig. 30. — SNR G41.1-0.3. Contours from Chandra observations: levels are 2, 4.3, 6.5, 8.8 and 11 X 10 ^ photons/cm^ /sec/pixel. The 
cross represents the location of the X-ray source J1907. 5-1-0708. Region used for partial photometry is also indicated in the figure (in white). 
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Fig. 32. — SNR G54.4-0.3. Contours from VGPS 21 cm observations: levels are from 9.5 to 14, in steps of 1.5 Jy/beam. Regions used 
for partial photometry are also indicated in the figure (in white). 
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Fig. 34.— SNR G304. 6+0.1. Contours from MOST observations: levels are 0.06 to 0.5, in steps of 0.11 Jy/beam. 
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Fig. 35.— SNR G310. 8-0.4. Contours from MOST observations: levels are 35, 58, 126, 240 and 400 mjy/beam. 




Fig. 37. — SNR G332.4-0.4. Contours from Chandra observations: levels are 0.8, 1.1, 2, 3.4 and 5.5 X 10"'' photons/cm^ /sec/pixel. The 
cross shows the location of the X-ray source 2E 1613.5-5053. Region used for partial photometry is also indicated in the figure (in white). 




Fig. 39.— SNR G337.2-0.7. Contours from Chandra observations: levels are 0.8, 1, 1.6, 2.6 and 4 X 10"'^ photons/cm^/sec/pixel. The 
cross shows the location of the X-ray source J163931. 4-475019. 




Fig. 41. — SNR G344. 7-0.1. Contours from MOST observations; levels are 16, 21, 35, 58 and 90 mjy/bcam. Circle represents the location 
of a catalogued IRAS point source corresponding to the diffuse emission peak. 




Fig. 43. — SNR G348. 5-0.0 (south-central countours) and the larger SNR G348. 5-1-0.1. Contours from MOST observations: levels are 
0.3, 0.7, 1.1, 1.6 and 2 Jy/beam. Diamonds represent maser locations. Regions used for partial photometry are also indicated in the figure 
(in black). 
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Fig. 45. — SNR G349.7-0.2. Contours from Chandra observations: levels are 0.3, 0.7, 1.9, 3.8, 6.5 and 10 X 10"^ photons/cm^/sec/pixel. 
Masers are represented by diamonds and the cross marks the location of X-ray source J171801. 0-372617. 
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